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TECHNICAL MEMORANDUM 1250 


SUBSONIC GAS FLOW PAST A WING PROFILE * 

By S. A. Christianovich and I. M. Yuriev 

The use of the linearized equations of Chaplygin to calculate 
the subsonic flow of a gas permits solving the problem of the flow 
about a wing profile for absence and presence of circulation 
(reference l). 

The solution is obtained in a practical convenient form that 
permits finding all the required magnitudes for the gas flow (lift, 
lift moment, velocity distribution over the profile, and critical 
Mach number) . This solution is not expressed in simple closed form; 
for a certain simplifying assumption, however, the equations of 
Chaplygin can be reduced to equations with constant coefficients, 
and solutions are obtained by using only the mathematical apparatus 
of the theory of functions of a complex variable. 

The method for simplifying the equations was pointed out by 
Chaplygin himself (reference 2). Tsien (reference 3) applied simi- 
lar equations to the solution of the flow problem and obtained a 
solution for the case of the absence of circulation. He did not 
succeed in obtaining the flow about the profile in the presence of 
circulation. 


1. EQUATIONS OF MOTION 


The equations of the two-dimensional adiabatic irrotational 
flow of a gas are of the form 

d(pu) d(py) _ n 

Sx Sy ' 


du 

Sy 



w 2 a 2 
~ Til 


a* 2 (X+l) 
2 (X-l) 


( 1 . 1 ) 


*" 0 btekanie Krylovogo Prof ilia Pri Dokriticheskoi Skorosti 
Pot oka." Prikladnaya Matematika I Mekhanika. Yol. 11, no. l, 1947, 
pp . 105-118 . 
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where w is the magnitude of the velocity, p the density, a 
the velocity of sound, a* the critical velocity, and u and v 
the projections of the velocity vector on the x- and y-axes. 

From these equations it follows that the expressions 


d,cp B ~ dx + ~ dy 
a* a* 


( 1 . 2 ) 


d\|/ = — H—X_dx + 5 !idy 

p 0 a* P 0 a * 

are total differentials; is the density of the gas at w=0. 

By solving these equations for dx and dy, there is obtained 


dx « °-° g -i dcp - fe * d\|/ 
X p X 

dy ss gjgJ. dcp+ gpjLi d\j/ 


(1.3) 


where X - w/a* . 

From the conditions expressing the fact that the right sides 
of equations (1.3) are total differentials 

c) sin _ S_ /Pq cos 

a \|/ * = a 

d cos 3 i Po sin 

2 ___ S^P ~ 

the equations of Chaplygin are obtained 

1 

aii if 

35 = * V 1 ~ xTTV 

ii , 1 ;?£ /, _ x-l ax 

acp “ * \ x+i / 
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ji 


By Introducing in place of the variable A the variable 
defined by the equation (reference 1) 


ds = 


1-A 


dA 


1-A (J0-D/(X+1) X 


(1.4) 


the equation of Chaplygin is transformed to the form 


ds 

15$ 



dd 

d^P 



(1.5) 



The Mach number is M = w/a. The constant of integration in 
equation (1.4) is chosen with the condition that e s /x = 1 for A-»0. 

The function e s is^ denoted in that which follows by X . 
Tables of the values of A and v/E as functions of X are given 

in reference 1. 

If in equations (1,5) (j. and y are taken as the unkown 

variables connected with s and d by the relation 

S - id a f ([i. + iv) 

where f (jjl + iv) is an analytic function of p. + iv, equations 
(1.5) can be replaced by the system, of equations' 

ds dd 

5T = ’as* 

( 1 . 6 ) 

dB dd 

Tu = 




and 
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& _ 

du 



(1.7) 


The coordinates x and y are determined from equations (1.3). 
The functions satisfying equations (1,6) and (1.7) convert the 
right sides of equations (1.3) into total differentials. 


2. APPROXIMATE EQUATIONS 

It is assumed as an approximation in equations (1.5) that 

VKs \/Kao where \/lC°o is the value of \/K corresponding to 
X = Xoo • 

Let 


\/ m \JT 

Equations (1.7) then assume the form 

Sep = M 

dp dv 

dep _ _ d^ 
dv dp 


( 2 . 1 ) 


Because of this simplification of equations (1.7), the right 
sides of equations (1.3) will no longer he total differentials 
after the solutions of the system of equations (1.6) and (2.1) 
are substituted in them. It is necessary to make changes in 
equations (1.3) to compensate for the simplification of equations 
(1.7). 

Then 


s* = X 2a -i- ; 

ds p \TZ P \/kZ 



Moreover, 


t 
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Whence 


k 

d / £o l\ _ _ \/~K _ a,/'Kos 

ds \o 7y X “ X 




1 

X 


i = + 0 2 e-B 


^0 1 
P x 


(C 2 e- s -C 1 e s ) \/K^ 


( 2 . 2 ) 


The values of the constants are chosen from the condition 
requiring the exact satisfying of equations (2.2) for X = Xqq. 


C 


1 


Xco i A _ i N 

2^x7 \ \/ l-M j) 




(2.3) 


Substituting in equations (1.3) the corresponding expressions 
from equations (2.2), the following equation is obtained: 

dx + idy - (dqvid^F) - C 2 e _s+i * (dcp+id\p) (2.4) 


The right side of equation (2.4) by virtue of equations (1.6) 
and (2.1) is a total differential. 

Consider the Jacobian D(x,y)/D(cp,i0 : Making use of equation 

(2.4), it is found that 

££4l = C p 2 ®- 20 - C n 2 e 2s 
D(cp,\lr) 2 X 



6 


NA.CA EM 1250 


In order that the determinant should not be zero, it is 
necessary that the following inequality be satisfied: 



This condition is a consequence of the modification made in 
the equations of motion; however, it is not essential. This 
inequality is generally satisfied if the fundamental inequality 

A*^ A is satisfied where A_„ = 0.7577 ... corresponds to A = 1. 
cr or 

3. FOEMUIATIOH 01 THE PROBLEM 

In the plane z - \i + iv, a closed contour about which is a 
flow of an incompressible fluid Is considered. The direction of 
the flow will be taken to coincide with the direction of the 
x-axis and the velocity at infinity A® • Let f = f (z) be the 

complex potential corresponding to this flow and set 


8 - U » log || (3.1) 

The solution of equations (1.6) as determined by equations 
(3.1) Is to be treated as the solution of the problem of the flow 
about a profile in the |i,u-plane of a fictitious stream of an 
incompressible fluid. The magnitude of the velocity at each 
point of this flow Is equal to A and the single of inclination 

of the velocity vector to the x-axis Is equal to dJ thus 


df 

dz 


= Ae-i* 


The complex potential may be represented in the form 


f 


( l * 


1 

i 




= x (I) 


z 


(3.2) 
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where z = X (I ) is a function mapping the outer region -with, respect 
to a circle of unit radius in tlie \ -plane on a region outside the 
contour considered. This function satisfies the conditions 

XH = 00 

t =1 

t =® 

The last condition determines the scale for the contour in 
the plane n + iv. 

The expansion of z = X (£ ) in the neighborhood of the point 
at infinity has the form. 



z = | + a„+^i + f^.+ . . . 

0 T 1 2 


(3.3) 


Hence in the neighborhood of the point at infinity 


dz 


8.T —1 


1 + Z i + -±. - 

_ 2iti T ^2 


+ • 


(3.4) 


The circulation of the velocity T ■ in the flow of an incom- 
pressible fluid is equal to f = 7%o- 


■ In the | -plane, the velocity potential cp and the stream 
function will be considered also as the velocity potential 
and the stream function of an incompressible fluid about a cir- 
cular cylinder. 

The complex potential 3? = cp + i\|/ is introduced and repre- 
sented in the form 

r -X, (i * i + & ij (3.5) 

where y is as yet an arbitrary parameter determining the cir- 
culation of this flow. 
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Tn the solution of the exact equations it is necessary, in 
place of the function (3.5), to find the exact corresponding solu- 
tion of the linear system (1.7). In reference 1, this solution 
is constructed by the method of successive approximations. 

Equation (2.4), in the notation assumed, may be written in the 

form 


+ ldy = 02 3J7S ® + °i f§ ® (3 - 6) 

By considering all magnitudes as functions of the variable \ 
the following expression is obtained: 

“ + ldy - °2 d * + c i if 1*1 < 3 - 7 > 

or transforming it again 

dx + idy = C 2 dz + G x ^ ^ (ff) dz (3.8) 

Thus, knowing the function z =X(!)> it Is possible to make 
use of equation (3.7) or (3.8) to determine x and y with the aid 
of quadratures. 

The conditions the functions z = X ( I) and equation (3.5) must 
satisfy in order that the obtained solution correspond to a phys- 
ically possible flow are considered. 

In order that the solution have a physical sense, it is nec- 
essary that the region in the x,y-plane corresponding to this ■ 
solution be a single-sheet region. 

A necessary condition for this condition is that the Inte- 
gral of dx + idy should be zero over any closed contour in the 
region j l| > 1» 

Because the expressions on the right of equation (3.7) are 
the Bum of analytic functions of | and | having no singularities 
in the region ] £| >1, the corresponding Integrals do not 

depend on the shape of the path of integration and are equal to 
zero in going around any closed contour not including the circle 

lei = i. 
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The values of the integrals along any closed contour, including 
the circle J|| = 1, are equal; consider the integral of dx + idy 
along a cirole of large radius with center at the origin of 
coordinates* 


In the neighborhood of the point at infinity 


1 

df/dz 





whence 


dF/d| 
df/dz = 


1 


2UL- 

2ni 



+ 





Moreover 


df cLF ~ 2 
dz d| = Ao ° 


2iti I 



By substituting in equation (3.7), integrating over the cir- 
cle jll = E, and passing to the limit as B-*» 

J dx + idy « C 2 (7 - 7) + (7 + J) 

Keplacing C 1 and C 2 by the corresponding values from 
equations (2.3) gives 


J dx + idy = 

Thus the condition that to any closed contour in the region 
||[ > 1 sho ul d correspond a closed contour in the x,y-plane con- 
sists in the fact that the previously undetermined value of 7 in 
the function (3.5) must be taken equal to 




(3.9) 
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This value of y is proportional to the potential in the 
passage around the closed contour in the a:,y-plane and therefore 
. is proportional to the circulation of the actual flow. 

■When equation (3.7) in the neighborhood of the point at infin- 
ity is considered, it is found that the ellipse 


x 



J - 7 0 


Aoo 




i + 


Mob 



(3.10) 


corresponds to the circle of radius B in the plane £ = p + ip. 

Thus the neighborhood of the point at infinity in the £-plane 
is uniquely mapped on the region of the point at infinity in the 
x,y-plane. 


4. INVESTIGATION OF CONTOUR IN X,Y-PLANE 


Equation (3.7) is considered near the oontour j|j « 1. 
derivative of the function (3.2) becomes zero on the circle 
at the two points 


e 


01 


•«0 


t 


02 


e i(* - e 0 ) 



The 

let = i 

(4.1) 


Then 


df ~ (I - Soi) (6 - e02) (4.2) 

d| = Aco |2 


In an analogous manner, the derivative of the function (3.5) 
where y is determined by equation (3.9) can be represented in 
the form. 
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where 


- eIS * 


1^2 = 


sin 8 . 


4« 




1 - ‘ 


Equation (3.7) will therefore he 

~ 2 , 


(4.4) 


dx + idy - C 2 fcfcjfcjgx- M At + C]> == d£ (4.5) 


In the absence of circulation, for example, for 7 <* 7 = 0, 


dx + idy =. C 2 X' (I) d | + CjXl dt 7 -J :) 2 . 


(4.6) 


It is shown hy equation (3.7) or (4.5) that the points of the 
contour corresponding to the zeros of the functions df/dz and 
dF/d| will, generally speaking, be singular points of the contour. 
The presence of these singularities is a consequenoe of the modi- 
fication of equations (1.7). 


In considering in greater detail what takes place in the 
neighborhood of these points, let df/dz in the neighborhood of 
| = | Q have an expansion in the form 

l 0 + l l " W + ’ * * _ 


df 


S- « -5o>‘ 


where 0 5 : k^$ 1. To a cusp of the profile in the z-plane (for 
example, the tail point of the profile with zero angle) corresponds 
the value k = 0 , and to a regular point of the contour corresponds 
the value k = 1. In the neighborhood of the point \ = | Q (from 
equation (3.7)) 


(I - ! 0 ) k 


d| + 


dx + idy 
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where A is a certain coefficient; thus, 
for k<l 


* - *0 - 1 (y - y 0 ) - rr^ (t - y 1 '* + • • • 

for k = 1 

* - + 1 (y - = A 1o s (£ " t 0 ) + * * * 

It follows that an angular point of the contour in the 
x,y-plane will correspond to the point f «= this angular 
point with angle jc(l-k) is turned concave to the flow; for 
k = 1 the vertex of the single approaches infinity. 

The point I = at which dp/d £ becomes zero is now con- 
sidered. Let df/dz at this point also become zero «-nri the expan- 
sion of this function into a series in the neighborhood of this 
point take the form 

df ft t \ l 

The point at which 1=0 is a regular point of the con- 
tour in the plane z = jx+iu. The point at which 1>1 >0 is an 
singular point of the contour. At this point, the angle is con- 
cave to the region outside the contour. 

In the neighborhood of the point 1 = 1*, 

dx + idy = B ( t - | m ) lm 1 d | + . . . 


“0 + ” l 0 ) + • • • 


and therefore 


x - x Q + i (y - y 0 ) = -JL_ ( £ - ^ ) 2 “ l + . . . 

To this point in the x,y-plane also corresponds a cusp. For 
1=0 this point is a cusp with its sharp point faoing the flow; 
for l >0 this point is angular with angle jt(2-1); for 1=1 
it is an ordinary point of the contour. 


1258 
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5. FLOW ABOUT AH ABC OF A CIRCLE 

Before proceeding to the consideration of the flow about a 
contour of arbitrary shape, the example of the flow about a con- 
tour corresponding to an arc of a cirole in the plane z * p. + iv 
should be considered. The function 


T T,2 

z * f + ki + ± Z £ . 

4 .fTH 

where k < 1 maps the region outside the "out" along the arc of 
the circle shown in figure 1 into the region outside the circle 

|l| * 1. 

The derivative 

iz = 1 _ 1 - k 2 

d? (I + ki) 2 

becomes zero at the points 

! 01 « - ki -V L - k 2 
l 02 * - ki a / 1 - k 2 

corresponding to the ends of the arc. The angle 0 Q is therefore 
determined by the equation 

sin 0 q *» - k 


Then 


dz 

n 


(I ~ Iqi) (I ~ I02) 

(| + ki) 2 


Only the case of flow for which the velocities at the ends 
of the arc are finite can be considered. In this case. 


df 7 (E+ki) 2 

d£** w |2 


2rti 


2ki 
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Thus 





By integrating this equation, 



1258 



MCA TM 1250 


15 


The constant of integration C is chosen so that for Moo-^0 
the initial arc of the circle is obtained. For Moo-^0, 
equation (5.1) becomes 

1 v 2 

x + iy = f + kiC + sLzJL. 

& | + ki 

It is assumed that C =5^/^ 

The ends of the contour obtained correspond to the points 
and anji i ' iie P oiXL,fcs I ol an ^ \ 02 become ordinary points 

of the contour. The tangents at the angular points* of the contour 
coincide at the upper and lower sides. With an accuracy up to 

k 2 for x and k^ for y, the coordinates of the end points of the 
arc will be 



The maximum thickness of the contour is equal to 



The arc corresponding to the case k * 0.04, M w = M cr = 0.803 
is shown in figure 2. 

If in the expression for x there are retained only the mag- 
nitudes of the first order and in the expression for y only the 
magnitudes of the order of k, the profile will be an arc of a 
parabola 


y 





■) 


s imi lar to the initial arc but with an aocuracy up to higher orders 
In figure 2 this parabola is shown by a dashed line. 
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Still another example of an irrotational flow about a contour 
corresponding to an ellipse in the plane z = m- + iu shall he 
considered . The function 



where 0 ^ r 1 
the ellipse 


maps the region 



on the region outside 



(1 + r 2 ) 2 


+ 



)1 - r 2 ) 2 


1 


After integration equation (3.7) gives 


x + ly . 0 2 (f f)+ o log + S + -||] 

Because (5.2) 

lim (l + r 2 (r 2 - 2) log i,zS, + £ + 1 V £ + i. - -i. 
r,0 ^ 2r3 ® £ * r £ 3| 3 


equation (5.2) for r = 0 assumes the form 


I + 


If € — 1-r Is denoted and in equation (5.2) magnitudes of the 
order are neglected 


x + iy = Or 


(t ♦ o*. 8 (1 ♦ 4^ 
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The profile represented "by this last equation is very close 
to the initial ellipse. 


6. FLOW ABOUT AN ARBITRARY CONTOUR 

In order to obtain in the z,y-plane a contour without angu- 
lar points, it is necessary in the auxiliary plane of the variable 
z = ij. + iv to be given a contour having angular points chosen in 
a special manner. 

A certain auxiliary, function is constructed that enables the 
obtainment of the necessary singularities on the contour in all 
cases considered. 

In the plane of the complex variable £ , a region outside 
the contour shown in figure 3 is considered. This contour con- 
sists of the circle D and the cut along the arc of the cir- 
cle D' tangent to D at point c. The points a and b are 
the angular points of the contour. The region considered is 
mapped the region outside the unit circle (fig. 4) so that the 
points a, b, and c correspond to the given points A, B, and C 
on this circle and the derivative of the mapping function at the 
point of Infinity is equal to unity. The radii of the ciroles D 
and D’ and the length of the arc abc are then determined. 

The required mapping function is determined by the following 
equations : 


T) _ _ i £ + 

\ " to 

u = tj + (n* - ti 0 ) log (n* - t0 

£ _ _ ± ^0 + 1 Sc u + 

+ i Uos u. - Uco 

where ■ 

= - i + (n* - n Q ) log Ol* + i) 
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where | Q , and | c are the coordinates of the points A, B, and 
C In the plane f ; and t| # and Hq are the coordinates of the corre- 
sponding points in the plane T) . 


For the derivative, 



1 + 


(8o-t«)(t.to) , 

(U - tcj (to - SoJ 


1 

l - So J 


-2 


( 6 . 1 ) 

The function considered is denoted by & s g(|). As is seen 
from, equation (6.1), this function satisfies the condition that 

g'(oo; = 1. 


The function g(£), as follows from its geometric sense, is 
regular for |£| > 1 and on the circle j|| = 1 has two singular 
points l = £q and £ = |*. 

At the point | = £ c , the derivative d£/d£ is finite 

The contour I in the p., t-plane is considered as having no 
angular points. Let the function z « X(l) map the region out- 
side the contour onto the region outside the circle 1 = 1. 
let £ and £ Q2 he points on the oirole 1 1] = 1 corresponding 

to the critical points for the flow considered about the contour, 
and £ -j. and te the corresponding critical points for the 

flow about the circle determined by the function F(£). By 

g 1 (t) and g (£), the values of the function g corresponding to 
2 

the points £ Q , £ #1 and £ Q2 , £ * 2 are denoted. The values £ cl 

and i C 2 may be chosen arbitrarily. Instead of the contour E in 

the n,u-plane, the contour E* is considered, vhich bounds the 
area on which the region |||> 1 is mapped by the function deter- 
mined by the equation 

*£ -x' (?) 8l ' (5) g 2 ’ ({) 


( 6 . 2 ) 
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To this contour in the x,y-plane will correspond a contour 
■without angular points. 


The functions g^’ and gg ' 


may he represented in the form 


Si 


01 


h 


1 


S 2 


t ~ t02 
t ~ 1*2 


h 


2 


where h-^ and 
respectively, 

for computing the coordinates in the x,y-plane, then assumes the 
form 


hg become zero at the points 
as tlog( | Equation (4.5), 


I = and 


e.t 


*2 

which serves 


A 


* 


dx + Id, = 0 2 x' (S) ds 0 X (••»> 


On integration, the singularities corresponding to h^ and hg 
vanish. 


In the x,y-plane a closed contour without angular points is 
obtained. The actual construction of the contour E* Is unnec- 
essary. For constructing the contour in the x,y-plane and com- 
puting the velocity field, only the derivative (6.2) is required. 

The contour in the x,y-plane will approximate the initial 
contour in the plane z = u + iu and will approach coincidence 
with the latter as M^O. The greatest difference in the con- 
tours will be near the critical points. By a rational choice 
of the magnitudes | ol and l c2 > the difference can be made a 

minimum. The oontour shown in figure 5 is that one corresponding 
to the case where the function z = X( £) is z = t where 


= 0.333 

7 - C = 1 
7 

e . = iso 0 

cl 

e c2 = 30° 
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It Is evident that the indicated method of constructing the 
function g is not the only one possible ; it was discussed only 
from considerations of simplicity. It is not difficult to define 
a function g making possible the construction in the x,y-plane 
of a contour with any given angular point. 


7. APPROXIMATE CONSTRUCTION OP CONTOUR 


It was previously shown how to obtain in the x,y-plane a 
contour without angular points. The effect the angular points 
of the contour will have, if they exist, on the flow is now con- 
sidered. For this purpose, in the £ -plane two auxiliary cir- 
cles are drawn with centers at the points £ 01 and | passing 

through the points ( ^ and | #2 * Tb® closed contour formed by 

the arcs of the circle 1 1| = 1 and the arcs of the auxiliary cir- 
cles are mapped into a closed contour in the x,y-plsne. 

By making use of equations (3.2) and (3.5) and noting that 
|f| = 1, it is seen that 


|| = 2i (sin 9 - sin 0 Q ) 

41 . || -t - |) 


= 2iA (sin 9 - sin e„) M = e -21 * 

cif “ * dz 


Then 


r 2 


dx+idy = -f + 4C n A 
L 2 sin0-sin© n 1 00 


il 

dz 


(sin0-sin0Q) (sinS-sine^Jj' dz 


(7.1) 


In particular in the absence of circulation 
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{' 


dx + idj = < Cp + 4C. Xqo 


dz 


sin 2 © 


} 


dz 


The coefficient of dz is nearly constant for points removed 
from | and | • The contour in the x,y-plane therefore approx- 
imates very closely the contour in the z-plane. 


The radius of the previously mentioned auxiliary circles is 
equal to 


p = 2 eln z 2 (vi - m j * ') * e o Qrhg ~ x ) 

On an auxiliary circle A | | | pe^ e ; therefore 

dA£ = ip ei £ de 

By expanding dz/df in a series in the neighborhood of the 
point i = it can he seen that 

§| ** a + bA | + . » . 


where 


a = 



e o + 



With an accuracy up to magnitudes- contai ning P t o^ the first 
degree 


dx + idy 



- ]^tg + ipe ie Cgade + 0 (p 2 ) 


(7.2) 


By maVir i g use of equations (7.1) and (7.2 ) } in the x,y-plane 
the closed contour is constructed. The part of the contour corre- 
sponding to the auxiliary circles is not a streamline. 

The details of the flow near the singular points may be seen 
in figure ^ 6 , which gives the mapping of the circle for = 0.333 
and Cy = 7 = 1. 
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The solution obtained, may he interpreted as the flow about a 
contour with suction of a small mass of gas through the surface of 
the contour near the critical points. It can he noted that the 
value of p is always very small; furthermore for very large 
values of the lift force Q ^ is small, for example, 0 Q Z 4.5° 

for C_ » 1. For large values of C the value of M is small 

7 7 09 

because must not exceed M c;j j for small values of C y the 

values of 0* are small. 

Taking into account the fact that the singular points are 
looated in the neighborhood of the critical points where the veloc- 
ities are near zero and that p is in all cases of very small mag- 
nitude, the practical result is obtained that it is possible in 
many oases to neglect the effect of the singular points on the flow 
about the remaining part of the contour and to carry out its con- 
struction without the use of the function g. 


8. THEOREM OF JOUKCWSKI, COMPUTATION OF THE MOMENT 
OF THE AERODYNAMIC FORCES 


In order to compute the forces exerted by the flow on the 
airfoil, use is made of the well-known formulas obtained in the 
application of the law of conservation of momentum to the mass 
of fluid contained between the contour of the wing and a certain 
dosed contour E which includes the wing contour. 


?x 

P* 


-/W -f, 

E n E n 
-J pdx -J ] 


pu (udy - vdx) 
> 

pv (udy - vdx) 


( 8 . 1 ) 


where P and P v are the projections on the x- and y- axes of the 

pressure forces of the flow on the wing. The integrals in equa- 
tions (8.1) are taken in the direction for which the passage about 
the oontour E is effected, that is, in the counterclockwise 
direction. 


By applying the law of conservation of momentum for com- 
puting the aerodynamic moment of the wing profile. 
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*J p (xdx + ydy) - J' P (xv - yu)(udy - vdx) (8.2) 
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If 


By making use of equations (3.7) and (6.2), the integrals 
entering equations (8.1) and (8.2) may he transformed into inte- 
grals taken over the corresponding contour in the £ -plane. 'The 
contour E is so chosen that to it in the | -plane vill corre- 
spond a circle with center at the origin. All the magnitudes 
entering under the integral sign in equations (8.1) and (8.2) are 
expanded into series in the neighborhood of the point at inf ini ty 
in the £ -plane. From equation (3.7), 


dx + idy » C 2 ■( 1 + JL- ( ■■ X, .. - 1 ) + 

2 1 2 *i£ \\A - M, 2 ' 



On integrating, this equation becomes 


+ 



x + iy 




4it 2 


f [■ 1 + ^'^VTT^] + - • -} + ao ' 


+ a Q i 


t " 

where Sq + Sq i is the constant of integration and a, = + 

i7^ is the coefficient of | ” 2 in the expansion of df/dz in 
the neighborhood of the point \ ® oo. 

Expanding u, v, p, p, x, y, dx, and dy in powers of l/r 
and integrating, the following expression is obtained: 
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(8.3) 


(a") 


E*y — ~ Pocr ' W C0 r 



24 


MCA 01 1250 


The theorem of Joukowski Is thus valid also for the flow 
about a body of a gas with large subsonic velocities. This 
theorem was proven earlier in reference 1 by making use of the 
exact equations. 

In a similar manner , the equation for the moment of the 
aerodynamic forces is 


M 


v 2 o, 2 / 

p °° « (2*7. 

Vi \ 1 




(8.4) 


The constant bq’ determines the position of the point 

relative to which the moment is computed. The moment for the 
corresponding profile in an incompressible flow is equal to 


M “ p oo V 2 ( 2 *7 ;l “ 7P(j) (*L = a **J 

where is the real part of the free term in the expansion of 

the function z «=X* (|) in the neighborhood of the point | =e» . 
Let 00* The c3lord of" the deformed profile in the com- 

pressible flow will be denoted by b and the chord of the pro- 
file in the corresponding incompressible flow by b°; in gen- 
eral the magnitudes referring to the incompressible flow will 
be denoted by a superscript circle. 

The moment coefficient for an incompressible flow C ° is 

m 

given by 

0 2jC7i “ 70 0 


The moment coefficient for a compressible flow C m is given 
by 


-m 
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C_° = yfb° , and therefore 

v 



^Cy° 

vO * 



The origin, of coordinates is chosen at the nose of the profile. 

Then 



is a function of the focal distance from the nose of the profile 
and the expression for the moment C^ 0 may he given in the form 


where 



o o 

ac/ s' 


+ CL 


m 



2 * 7 -, 

‘1® 


is the moment of the profile in the incompressible flow relative 
to the focal point. 


The final expression for C m can be given in the following 

■ form: 


m 



+ C, 


m 


’] 


a 

OB 
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1 - b 2 


(8.5) 


Translated by S. Eeiss 
National Advisory Committee 
for Aeronautics. 
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